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A spectroscopic survey on the multiplicity of high-mass 
stars 
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ABSTRACT 

The formation of stars above about twenty solar masses and their apparently high 
multiplicity remain heavily debated subjects in astrophysics. We have performed a 
vast high-resolution radial velocity spectroscopic survey of about 250 O- and 540 
B-type stars in the southern Milky Way which indicates that the majority of stars 
(> 82%) with masses above 16 solar masses form close binary systems while this 
fraction rapidly drops to 20% for stars of 3 solar masses. The binary fractions of O- 
type stars among different environment classes are: clusters (72 ± 13%), associations 
(73 ± 8%), field (43 ± 13%), and runaways (69 ± 11%). The high frequency of close 
pairs with components of similar mass argues in favour of a multiplicity originating 
from the formation process rather than from a tidal capture in a dense cluster. The 
high binary frequency of runaway O stars that we found in our survey (69% compared 
to 19 — 26% in previous surveys) points to the importance of ejection from young star 
clusters and thus supports the competitive accretion scenario. 

Key words: stars: formation - stars: early-type - binaries: spectroscopic - techniques: 
spectroscopic 



1 INTRODUCTION 

The most massive stars, historically classified as O-type 
stars, are rare objects with masses above 16 Mq, typically 
found at large distances from the sun; early B stars 
(M > 8Mq) also contribute to the group of high- mass 
objects. The formation of these O- and B-type stars 
is still under debate and can b e explained by various 
models (|Zinnecker fc Yorkd 1200 it). Both observations of 
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formation scenarios, particularly if the multiplicity fraction 
were a function of stellar mass. While many massive stars are 
found to be part of binary or multiple systems, comprehen- 
sive statistics on close binaries is still missing. The smallest 
separations are expected to be around 0.2 AU, resulting in 
orbital periods of a few days only. Below this minimum dis- 
tance the binary system will merge to form a single object. 

The vast parameter space in possible orbital periods and 
mass ratios requires different, partly overlapping comple- 
mentary methods th at have their own lim itations and their 
observational biases (ISana fc EvandlioTll ): High-resolution 
imaging like adaptive optics or interferometric techniques 
serve for systems with wider separations and mass ratios 
q = M2/M1 between 0.01 and 1 while high-resolution spec- 
troscopy is biased toward finding close companions and 
those that are a significant fraction of the primary's mass 
(q > 0.1). We note that the inclination of the orbit with re- 
spect to the observer and the eccentricity are other crucial 
limitations for the spectroscopic detection of close multiple 
systems. 

An adaptive optics survey of about a third of the 
known gal actic O stars revealed visual companions in 27% of 
the cases (|Sana fc Evansll201ll ). Speckle interferometry for 
of stellar systems is a crucial constraint for the various star most of the galactic O stars showed companions for 11% of 
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iKuiper et alj |2010| ) seem to favour the accretion sce- 
nario. Iloweve£ 2 _the high multiplicity amon g high- mass 
stars jPreibisch et al J 1 19991 : iMason et afll2009h might alter- 
nativ ely support a merging process o f intermediate-mass 
stars (|Bonnell. Bate fc Zinneckerlll99Sh . 

Massive binary stars are believed to be the progen- 
itors of a variety of astrophysica l phenomena, e.g. short 
gamma-ray bursts dEichler et alJ 1 19891 : iPaczvnskil Il99ll: 
Naravan. Paczvnski fc Piranl 1992), X-ray binaries ( Moff atl 
20081). millisecond pul sar systems and double neutron stars 
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the sample l|Mason et alj|200sh . Moreover, the same study 
claimed that 51% of the O-type objects are spectroscopic 
binaries (SBs), based on an extensive review of the liter- 
ature. This is in accord with a recent spectroscopic survey 
finding that among ~ 240 southern galactic O and WN stars 
more than one hundred star s show radial veloci ty (RV) vari- 
ations larger than lOkm/s (|Barba et al.ll201oT ). Recently, a 
high-resolution imaging campaign of 138 fields containing 
at least one high-mass star yielded a multiplicity fraction 
close to 50% |Mafz-A pcllaniz 20icf). In summary, the spec- 
troscopic binary frequency of high-mass stars so far observed 
and reported in the literature is moderately high, while the 
visual binary fraction is low. Little or no discussion as to why 
the multiplicity for O-type stars is high has been offered up 
to now nor have conclusions been drawn so far about which 
of the competing high-mass star formation models best ex- 
plains the hitherto known trends in the stars' multiplicity. 

There have been several surveys of B star duplicity in 
the past. In a sample of 109 B2 - B5 stars there were 32 
(29%) spectroscopic and 49 (45 %) visual binaries yielding a 
total binary frequency of 74% i|Abt. Gomez fc Levil fl990). 
A spectroscopic survey of 83 late B-type stars revealed that 
24% of the stars had companions with mass r atios greate r 
than 0.1 and orbital periods less than 100 days l|Wolfllll97St ). 
A speckle interferometry survey of the Bright Star Catalogue 
resolved 34 of 245 B stars in to binaries correspond i ng to a 
multiplicity fraction of 14% (McAlistcr et al ] |1987| , |1993D . 

Another speckle interferometry survey o f 48 Be stars 
revea led a similar binary fraction of 10% ± 4% l|Mason et al.l 
Il997h . A further comparison study, based on adaptive optics 
IR imaging and probing separations from 20 to 1000 AU for 
40 B and 39 Be stars derived the same qualitative result, i.e. 
that the multiplicity of B and Be stars are identical. This 
time, however, the binary fractions were 29% ± 8% for the 
B sta rs and 30% ± 8% for the Be stars (jOudmaiier fc Parrl 
2010). Finally an adaptive optics photometry and astrom- 
etry su rvey of 70 B stars revealed 16 resolved co mpanions 
(23%) (|Roberts. Turner fc Ten Brummelaaill2007h . In sum- 
mary, the overall multiplicity fraction of high-mass stars 
seems to decrease with stellar mass. 

In the present paper we investigate the multiplicity frac- 
tion in the stellar mass range of about 3 — 80 Mq and for 
mass ratios q > 0.2. 



2 OBSERVATIONS 

We have performed a comprehensive spectroscopic survey 
on a large representative sample of 249 O- and 581 B-type 
stars. The O- stars were taken from the Galactic O-Star Cat- 
alogue V.2.0 (|Sota et alj|2008l ) (GOSC) which is assumed to 
be complete to V ^ 8 mag. The spectral classification in the 
GOSC is curr ently under revi sion: A new version has been 
presented by (|Sota et al.ll201lT ) for 184 stars with B < 8 and 
north of 8 — —20°; a second revision is in preparation by the 
same authors. The B stars were selected from the HIPPAR- 
COS archive in a way that each subclass (B0 - B9) contains 
roughly an equal amount of stars; 50% of the B stars form 
a volume- limited sample with d < 125 pc. The distribution 
of visual magnitudes is displayed in Fig. Q] showing that the 
B-star sample contains on average brighter stars than the 
O-star sample. 




9-10 10-11 11-12 
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Figure 1. Distribution of V magnitudes for the samples of O- 
type stars (grey) and B-type stars (black) columns. The B stars 
are on average brighter which results in better S/N spectra and 
thus increases the chance for detecting fainter companions. 



Using the high -resolution spectrograph BESO 
l|Fuhrmann et all 1201 lh at the Hexapod Telescope at 
the Universitatssternwarte Bochum near Cerro Armazones 
in Chile we obtained 3632 multi-epoch optical spectra 
were obtained with BESO. The observing period started in 
January 2009 and is still going on. The spectra comprise a 
wavelength range from 3620 to 8530 A and provide a mean 
spectral resolution of R = 50.000. The entrance aperture 
of the star fibre is 3'.' 4. The integration time per spectrum 
was adapted to the published visual brightness of each star. 
It was our primary goal to monitor a large number of stars 
rather than to obtain a very high S/N for individual stars. 
As a consequence of this strategy potential companions 
fainter than AV ~ 2 mag are barely visible in our spectra, 
thus decreasing the chance for detecting SB2s. Converting 
this brightness difference of 2 mag into a mass difference 
we are sensitive to mass ratios q > (0.18 — 0.40) for 05 - 
09 stars and q > (0.43 — 0.55) for B stars. In other words, 
the detectable companions of an 05 star range from 05 to 
about B2, those of a B9 star from B9 to about A7. 

Additional spectra were collected during various observ- 
ing runs with FEROS at the ESO 1.5 m and the MPG/ESO 
2.2 m telescopes, both located on La Silla, Chile. Because 
BESO is a clone of FEROS the instrumental parameters are 
identical to those described above. These spectra cover a 
time span between 2006 and 2008. 

Finally we extracted 1009 O-star spectra and 133 B-star 
spectra from the ESO archive covering an observing period 
from 1999 to May 2010. These spectra were also obtained 
with FEROS and processed in the same way as the BESO 
spectra. 

Until today there are 2565 O- and 2209 B-star spec- 
tra in our archive. On average we collected ten spectra for 
each O-star with a minimum of two spectra for those ob- 
jects which were known to be multiple systems before and a 
maximum of 30 spectra in such cases where multiplicity was 
not immediately obvious. The spectra are separated by days, 
weeks and months. So far we have observed 540 stars from 
our B-type sample with an average number of four spectra 
per star. We stopped temporarily collecting data for those 
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Figure 2. Multi-epoch BESO spectra obtained between 
JD 2454878 and 2455295 for the double-lined O star binary 
HD 104649. The selected spectral regions show several Hel linos, 
Ha and the interstellar NaD lines. The change from a single- to 
a double-lined structure is clearly visible. 



Table 1. Spectroscopic results for 248 O- and 581 B-typc stars. 



Name 


Sp.T. 


No. 


Mult. 


HD 000560 


B9V 


10 


C 


HD 014951 


B7IV 


10 


SB2 


HD 017543 


B2IV 


10 


SB1 


HD 035588 


B2.5V 


6 


U 



B-stars where the multiplicity became obvious during the 
first two spectra. 

All data were reduced with a pipeline based on the MI- 
DAS package developed for the ESO FEROS spectrograph. 
A quantitative analysis was done with standard IRAF line- 
fit routines allowing the detection of line shifts in single-line 
spectroscopic binaries (SB1) or line deformation and/or sep- 
arations in double-line spectroscopic binaries (SB2). For the 
O- and early B-stars we used exclusively He lines for iden- 
tification; Hel (A5875 A) was particularly useful due to the 
nearby interstellar Nal doublet (AA5890, 5896 A) allowing 
for a sensitive verification of any relative line shift. For later 
B-types we had to rely primarily on hydrogen lines. A typi- 
cal set of multi-epoch spectra is shown in Fig. [2] 



3 RESULTS 

The data of our survey will be published in electronic for- 
mat. Table [1] shows the corresponding entries: Name, Sp.T. 
= spectral type, No. = number of spectra, and Mult. = mul- 
tiplicity status (C = constant RV, SB1, SB2, U = unknown) 
and some example results. 

We detect RV variations for about 82% of the O stars 
brighter than V = 8 mag while for the B stars the vari- 
ability fraction decreases from 61% (B0) to ~ 15% (B9) 
(Fig. [3}. For fainter O stars the cumulative percentage of 
RV variations diminishes to 71% (V < 9 mag) and 68% 
(V < 10 mag) reflecting the decreasing instrumental sensi- 
tivity. In general, our observed multiplicities are lower limits 
due to the unknown orbit inclination and the constraints for 
the detectable minimum mass ratio. 

The overall percentage of SBs for O stars is higher 




03.5-6 07-8 09 BO-1.5 B2-3 B4-6 B7-8 



Figure 3. Multiplicity vs. spectral type of O- and B-type stars. 
The black columns denote the total percentage of spectroscopic 
binaries, the grey columns display the percentage of double-lined 
spectroscopic binaries. The number N on top of the columns de- 
note the number of objects in each spectral bin. 



than found before in similar investigations and is due to the 
numerous spectra obtained for a single star, which reveal 
mainly variations within days and weeks that were not ob- 
vious in previous studies. As mentioned above [Mason et all 
(2009) report a general spectroscopic multiplicity fraction 
of 51% for the population of Galactic O stars. Individual 
nearby clusters were found to have binary fractions be- 
tween 0% and 63% with an average value of 44% (see e.g. 
ISana fc Evan! (120111) for an overview). 

Inspecting those 60 O stars that did not show any RV 
variations throughout the last six years, we find that ad- 
ditional 13 out of 29 stars observed t hrough adaptive op- 
tics mea surements |Mason et al.| [l997') or speckle interfer- 
ometry dMcAlister et alj 1 1993ft ' possess visual companions. 
These complementary data increase the total percentage of 
multiple systems for stars brighter than V = 8 mag to 91% 
leading us to suggest that basically all O-stars are mem- 
bers of multiple syste ms. This finding also g oes beyond the 
results summarized by ISana fc Evansl l|2011 ) who obtain a 
total minimum multiplicity fraction close to 70%. 

Among the SBs in the present study 82% of the spectra 
for O stars with V ^ 10 mag contain more or less separated 
multiple lines (SB2) reflecting that the majority of systems 
contains pairs of similar mass . This is in agreement with 
results bv lKobulnickv fc Frverl l|2007ft who found that mas- 
sive stars preferentially have massive companions. From our 
current time coverage, however, it is not yet possible to de- 
rive orbital periods and constrain semi-amplitudes of radial 
velocities and hence accurate binary component mass ratios. 

Due to the limited spectroscopic binary studies for B 
stars a comparis on between our work a nd pr evious investiga- 
tions is difficult. iKouwenhoven et al.l l|2005ft studied the bi- 
nary population in the Scorpius OB2 association. Summariz- 
ing the results from all techniques they find multiplicity frac- 
tions of - 80% and ~ 50% for B0 - B3 and B4 - B9, respec- 
tively. While the numbers for the early B-types are compat- 
ible with our study the multiplicity rate for the late B-types 
is slightly higher in Sco OB2; the latter difference might be 
due to the different observi ng techniques, because wi thout 
the adaptive optics data of IKouwenhoven et alj l|2005ft the 
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multip licity rate would drop by a bout 10% for the late B- 
types. iKouwenhoven et al.l (|2007l ) repeated the analysis of 
the primordial binary population in Sco OB 2: one set of 
spectr oscopic data in this study comes from iLevato et al.l 
(1 19871 ) who found that the binary fraction is at least 30% 
for all early-type stars but might be as high as 74% if all 
reported RV va riations were due to binaries . Another spec- 
troscopic s tudy jBrown fc Verschuerenll997h included in the 
analysis bv lKouwenhoven et al.l 1 2007h yields a similar range 
between 28% and 76%. However, the dependence on stellar 

mass has not been addressed in these stud ies. 

IWheelwright, Qudmaiier fe Goodwin! l|201Cf ) investi- 
gated the binarity of 25 Herbig Be stars with spectro- 
astrometry and derive a high binary fraction of 74%. 
This result, however, is in contrast to the AO study of 
lOudmaiier fc Parr! [|2010u who found a binary fraction of 
30% ± 8% for 39 Be stars. The same authors claim that the 
binarity of normal B stars is 29% ± 8% and thus identical 
to those of Be stars. In summary, the multiplicity of B stars 
is still a subject that needs further investigation; the exist- 
ing results - including ours - likely suffer from various bi- 
ases. We actually might expect to observe more binary stars 
amongst late type B stars compared to O stars, however, 
then their mass ratios must be larger than what is covered 
by our study. 

One may speculate whether the decrease of the bina- 
rity fraction between O- and late B-type stars is due to 
an evolutionary effect: an O-type star - whether on the 
main sequence or not - is much younger than a late B-type 
star on the main sequence and is likely more evolved due 
to its rapid hydrogen burning. If evolution "destroys" bi- 
narity then evolved stars (i.e. luminosity classes III and I) 
should have lower binarity fractions than stars on the main 
sequence. Fig. U shows the distribution of luminosity classes 
I, III, and V for the different spectral type bins in our sam- 
ple. Obviously those spectral types with the highest binarity 
fraction (03 - Bl) have the highest fraction of evolved stars 
while late B-Type stars with the lowest binarity fraction 
comprise only few evolved stars. This excludes the possibil- 
ity that evolution influences binarity. On the other hand, 
one can exclude a pure age effect as SB2 binaries are all 
very tight and hard to break up. They could dynamically 
interact in dense cluster cores and get ejected, but breaking 
up such "hard" SB2s is not very likely. 

The multiplicity of high-mass stars seems to depend on 
their environment. It was found that the binary frequency 
among O stars in clusters and associations is much higher 
than among field stars (which have no apparent nearby clus- 
ter) and runaway stars (O stars with peculiar RVs in excess 
of 40km/s or remote from the Galactic plane). The spec- 
troscopic binary fraction in clusters and associations, in the 
fiel d, and among runaways obtai ned in two different stud- 
ies (|Giesll 19871 ; TlVIason et al.lll998T l were 55 (61)%, 45 (50)% 
and 19 (26)%, respectively. The O stars in the present sam- 
ple can be assigned to clusters (44), associations (123), field 
(25), and runaways (54). For th e classification w e have used 
the designations in the GOSC |Sota et al.ll2008h refined by 
a recent investigation which turned 73 out of 93 "field stars" 
from the GOSC into runaway s originating from open clus- 
ters (ISchilbach fc Ro scr 2008). In those cases where stars in 
clusters (4), associations (3) or in the field (47) were classi- 
fied as runaway stars we have adopted this new designation. 




03.5-6 07-8 09 BO-1.5 B2-3 B4-6 B7-8 B9 

Figure 4. Luminosity classes for each spectral bin. 



The fractions of SBs within the individual groups are: clus- 
ters (72 ±13%), associations (73 ±8%), field (43 ±13%), and 
runaways (69 ± 11%); the errors are based on the number 
statistics. While our results for the clusters/assoc iations and 
the field are comp atible with previous studies l|Gieslll987l : 
iMason et al.ll 19981 ) our binary fraction among t he ru n aways 
(69 ± 11%) is significantly h i gher compared to iGiej (ll987T ) 
(19 ± 5%) and lMason etall (|l998l ) (26 ± 5%), respectively. 
If field stars are dynamically ejected from a cluster, one ex- 
pects an equal fra ction of singles and binaries expelled in op- 
posite directions IIHoogerwerf. de Bruiine fc de Zeeuwll2000l . 
l200ll : lGualandris. Portegies Zwart fc Eggletonll2004l ). 



4 DISCUSSION 

While the direct observation of disks around the earliest 
O stars will remain a fundamental challenge to test model 
computations for individual high-mass systems, binary star 
statistics offers other constraints on hydrodynamical sim- 
ulations of star forming clusters. Early calculations of the 
fragmentation of an isothermally collapsing cloud have al- 
ready shown that binary systems and hierarchi cal multiple 
systems are frequently obtained jLarsonl 1 19781 ). Although 
the variety of processes and th eir sequences are manyfold 
|Bate. Bonnell, fc Brommll2002l ) there is consensus that the 
accretion processes favor the formation of binaries with mass 
ratio s q ~ 1 for tho se systems with separations below 10 AU, 
(e.g. IClarkd (|2007h and references therein). 

Nevertheless, the simulations almost never produce 
binaries with q < 0.5. This appears to be a general 
problem of turbulent fragmentation calculations which 
seems to origin from gas accretion onto a p roto-binary 
jGoodwin. Whitworth fc Ward- Thompson! 120041 ). Indepen- 
dent of how the binary system was formed, the infalling 
material has a high specific angular momentum compared 
to the bina ry and thus will prefe re ntially accre te onto the 
secondary (|Bate fc Bonnelll l2004h . iBatel (|2000l ) even pre- 
dicts that binaries will evolve rapidly towards q ~ 1 re- 
gardless of their initial mass ratio. For the special case of 
massive stars recent theoretical calculations claim that close 
high-mass stellar twins can in principle be formed via frag- 
mentation of a disk around a massive protostar and subse- 
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quent mass transfer m such close, rapidly accreting oversized 
proto-binaries jKrumholz fc Thompsorj2007l ). This scenario 
provides a natural explanation for the numerous high-mass 
spectroscopic binaries with mass ratios close to unity. 

The high binary frequency for clusters found in the 
present study exclude random pairing from a classical IMF 
as a process to describe the similar-mass companions in mas- 
sive binaries. Our results make binary-binary interactions 
inside clusters very probable and thus can explain the high 
binary fraction among runaway stars. Because ejec tion re- 
quires a cluster origin of the binaries (|Kroupal l2001) it sup- 
ports the model of competitive accretion within the cluster 
environment. Finally, the small number of field O stars - 
only six "certified" stars remain currently in this category - 
suggests that probably all O stars are born in clusters. 

Radiation- hydrodynamic simulations show that, during 
the collapse of a massive prestellar core, gravitational and 
Rayleigh- Taylor instabilities channel gas onto the star sys- 
tem through non-axisymmetric disks. Gravitational insta- 
bilities lead to a fragmentation of the disk around the pri- 
mary star and form a massive companion. Radiation pres- 
sure does not limit stellar masses, but the instabilities that 
allow accretion t o continue lead to small multiple systems 
(|Krumholz et alj|2009h . 

The current study was a pure discovery project that 
aimed exclusively at the multiplicity statistics in the mass 
range 3 < M[Mq] ~ 80. In a next step we will study the or- 
bital properties and the spectral types of the individual com- 
ponents to obtain a statistically relevant archive for high- 
mass binary systems. Likewise, we are currently performing 
a photometric monitoring program to search for eclipsing 
binaries in order to determine accurate masses for O stars. 
From the fact that most O-type stars occur as similar mass 
binaries we conclude that their absolute magnitude calibra- 
tion has been overestimated; in the worst case the tabulated 
O star magnitudes are wrong by 0.75 mag at each wave- 
length. 
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